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ABSTRACT 


The paper presents and evaluates a force-derivative method which 
produces higher-order modal solutions to transient problems. These 
higher-order solutions converge to an accurate response using fewer 
degrees-of-freedom (eigenmodes) than lower-order methods such as the 
mode-displacement or mode-acceleration methods. Results are presented 
for non -proportionally damped structural problems as well as thermal 
problems modeled by finite elements. 


SYMBOLS 


C damping matrix 

E modulus of elasticity 

e spatial error norm (eq. (34)) 

K stiffness matrix 

K generalized stiffness matrix, conductance matrix 

M mass matrix, moment 

m subset of total number of degrees-of-freedom 

M generalized mass matrix, capacitance matrix 

n total number of degrees-of-freedom 

Q force vector 

Q generalized force or thermal load vector 

q modal coordinates of the second-order system 
t time 

u displacement 

x coordinate direction 

Y generalized displacement or temperature vector 
Z modal coordinates of the first-order system 
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Greek 


oq ith damped eigenvalue (eqs. (3) and (4)) 

[a] matrix of damped eigenvalues 

<I>i ith damped eigenvector 

[*] matrix of damped eigenvectors 

. (j>i ith normal eigenvector (eq. (21)) 
e time-integrated error norm defined in eq. (33) 

£2 matrix of frequencies squared (eq. (23)) 
to, ith circular natural frequency 

(o d j ith circular frequency of the damped free vibration 
A matrix of damping coefficients (eq. (23)) 
p density 

x dummy variable of integration, temporal integral limit 
(eq. (33)) 

Ci ith modal viscous damping factor 

Subscripts 

o initial 
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Superscripts 

a approximate 

(i) ith derivative with respect to time 

A matrix of reduced number of eigenvectors or eigenvalues 

or reduced number of modal coordinates 
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INTRODUCTION 


Transient thermal and structural analyses of complicated engineering 
problems which are modeled as discrete multidegree-of-freedom systems 
often require the solution of very large systems of equations. Reducing the 
order of such systems is highly desirable from the standpoint of increased 
computational efficiency. Some of the many methods for reducing the 
order of discrete multidegree-of-freedom structural dynamic systems 
include mass condensation methods (e.g., refs. 1 and 2) and reduced basis 
methods (e.g., refs. 3-7). The reduced basis methods use either a truncated 
set of basis vectors (e.g. eigenmodes, Ritz vectors, or Lanczos vectors) or a 
combination of basis vectors (e.g., eigenmodes and Ritz vectors (ref. 8)). 
Reduced basis methods have also been applied in solving transient thermal 
problems (refs. 9-11). Some problems, such as the thermal/structural 
analysis of space transportation systems or large space structures, require 
a large number of basis vectors to accurately represent the transient 
response. In addition, when singularities occur in the loading, convergence 
of a solution is not guaranteed. For most transient thermal problems, it 
takes a large number of eigenmodes to accurately model the sharp thermal 
gradients within the domain (ref. 10). 

When a reduced basis method uses the eigenmodes of a system of 
equations the method is referred to as a modal method. Two of the most 
widely used modal methods for transient structural analysis are the mode- 
displacement method (MDM) and the mode-acceleration method (MAM). 

It was shown in reference 12, that the MAM can be considered a higher- 
order modal method than the MDM and that the MAM converges to an 
accurate solution using fewer modes than the MDM (refs. 6 and 12). A 
method for generating improved or higher-order modal methods was 
developed in reference 7. This method was generalized to proportionally- 
damped structures (the damping matrix is a linear function of the stiffness 
and mass matrices) and evaluated for a variety of loads and damping 
ratios (ref. 12). This new method successively integrates-by-parts the 
convolution integral form of the solution and is called the force-derivative 
method (FDM) because it produces terms which are related to the forcing 
function and its time derivatives. The FDM was found to be more accurate 
than either the MDM or MAM. In particular, for problems in which there 
are a large number of closely-spaced frequencies (e.g., large truss-type 
space structures or multispan beams) the FDM is very effective in 
representing the important, but otherwise neglected, higher modes. Recent 
work indicates that the FDM, which was based on reference 7, has several 
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variants (refs. 13 and 14) for solving transient structural problems. One 
such variant is called the dynamic-correction method (DCM) (ref. 13) which 
is useful when a particular solution exists for a given forcing function. 

The purpose of the present paper is to extend the work of reference 12 to 
include non-proportionally damped structural systems and to evaluate the 
usefulness of the FDM in solving non-proportionally damped structural 
problems as well as thermal problems. Modal methods evaluated include 
the MDM, MAM, FDM, and the DCM. 
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THEORY 


First-Order or Damped-Mode Formulation 


The equations of motion, in matrix form, of an n degrees of freedom 
system, together with the initial conditions are given by 


Mii + Cu + K u = Q(t) 


( 1 ) 


u(0) = u 0> u(0) = u 0 


where M, C, and K are the mass, damping, and stiffness matrices of the 
system; u and Q are the displacement and load vectors, respectively, and a 
dot denotes differentiation with respect to time. 

Transforming eq. (1) to first-order form results in the following system of 
equations: 


MY+K Y= Q 


Y(0) = Y 0 (2) 


where 



and 
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, Y 
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Assuming a solution to the homogenuous form of eq. (2) as: 

a r t 

Y(t)= e 4> r (3) 

results in an eigenvalue problem 

a r M O r + K d> r = 0. (4) 

For distinct eigenvalues (a r ), the eigenvectors (<I> r ) are normalized such 
that 

T _ 

O r M <P r = 1.0 

and then (5) 

T _ 

<l> r k <t> r = -a r 

Equation (5) can also be written in matrix form as 


[<D] T M [<D] = [I] and [0] T K [<D] = [a] 


where [<t>] is a 2n by 2n modal matrix with its ith column equal to <b r an( j 
[a] is a diagonal matrix consisting of the a r ’s. 
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A solution to eq. (2) is assumed in the form of the following modal 
summation 


2n 

Y(t) = ]TO r Z r (t) 

r=l 


( 6 ) 


T 

Substituting eq. (6) into eq. (2) with premultiplication by O r results in the 
following, uncoupled, system of equations: 


T _ 

Z r -a r Z r = <D r Q 


Z r (0) = Z ro = 0 T r M Y 0 


(7) 


The solution to eq. (7) is 


a r t f CCr(t-T) T _ 

Z r (t) = Z ro e + J e O r Q(x)dx ( 8 ) 

0 


Hence, the solution of eq. (2) becomes 


Y(0 = 



r= 1 


a r t 

e 


f a 

+ Je 
o 


r(t-X) 


F r (x)dx 


where 


(9) 


F r (x) = <D T r Q(x) 
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If the forcing function has continuous derivatives, the convolution integral 
of eq. (9) can be integrated by parts to produce higher-order modal 
methods (ref. 12). For example, if it is integrated by parts once, the 
following expression results 


Y(t) = 



F r (0)l <M O r F r (t) 
ct r . a r 


r=l 


t 


+ 


Or J^t-T) 
a r 0 


F r (X)dT) 


( 10 ) 


If all the modes are used in the second-to-last term in eq. (10), this term 
can be written as 


2n i 

X-* r fe] Fr(t) = •t't>][^]['i>] T Q(0 =k‘q(0 (11) 

r=l 


if equation (11) is substituted into equation (10), the resulting expression 
is analogous to the MAM (ref. 6) and can be written as 
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_-l_ 

+ K Q(t) 


( 12 ) 


If the forcing function has continuous derivatives up to order N (C N ), the 
convolution integral of eq. (9) can be integrated-by-parts N-times, 
resulting in the following expression 



where the superscript (i-1) denotes the (i-l)th derivative. 

If all the modes are used in the last N-terms of eq. (13), they can be 
represented as functions of M and K (ref. 12) as follows 




Alternate Damped-Mode Formulation 


Equations 8 to 14 are analogous to those presented in reference 12, 
however, the present expressions solve a first-order system of equations, 
using the damped modes O r to decouple a non-proportionally damped 
system (the damping matrix is not a linear function of the mass and 
stiffness matrices). This means the order of the system of equations is 
doubled (2n). Equation 14 represents a means for developing higher-order 
modal methods than either the MDM or MAM, as demonstrated in 
reference 12, and is called the force-derivative method (FDM) because it 
produces terms which are related to the forcing function and its time 
derivatives. 

Equation 2 can also be considered to represent a heat conduction problem 

where, for that problem, M would represent the capacitance matrix, K 

would represent the conductance matrix, Q represents the thermal load 
vector, and Y is the vector of nodal temperatures. 

The MDM uses a subset, m (m < 2n), of the eigenmodes to reduce the size of 
the problem and solves for Zr using eq. (8) and substitutes these values 
into a reduced modal summation in eq. (6) to approximate the response, 
Y(t). The MDM can be classified as a zeroth-order method because it is 
equivalent to using the FDM (eq. (14)) with N = 0. An analogous form of the 
MAM uses eq. (12), and a reduced modal summation to approximate Y(t) 
and can be classified as a first-order method (N = 1 in eq. (14)). The FDM 
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uses eq. (14) with N>1. Reference 12 showed that the expressions obtained 
with four integrations by parts offer improved approximations to the 
higher, neglected, modes for several structural problems. 

The FDM (eq. (14)) can be derived using an approach similar to that used 
in reference 13 which results in a form which is more suitable for inclusion 
into existing thermal and structural analysis codes. A numerical approach 
can be derived, similar to that presented in reference 13, which 
approximates the forcing function as a piecewise differentiable polynomial 
and which numerically integrates the reduced system of equations 
(eq. (8)). For example, assume the forcing function is differentiable. 
Equation 7 could be differentiated twice to produce the following equations 


Ti 

Z r - (XrZr =<I> r Q 


(15a) 


(3) .. T - 

Z r - a r Zr =O r Q (15b) 


Re-arranging eq. (7) and substituting for Z r and Z r from eqs. (15a and 15b) 
results in 


z r(0 = - — <*>r Q(0 


^2<E> T r Q(t)- -^-o’rQO) 

ccf at 



(3) 

Z r (t) 


(16) 


Using eq. (6), the response can be written as 


1 3 


-1 ~ 

Y(t) = K Q(t)-KMK Q(t)+KMK MK Q(t) 


2 n 



r=l 


(17) 


The last term can be evaluated using eq. (8) and Leibnitz’s rule for 
differentiation of an integral to produce the following 


(3) 2 T — T - T - 

Z r (t) = cc r <D r Q(t) + a r O r Q(t) + 0> r Q(t) 

t 

a a r t 3 f a r (t-X) T _ 

+ Or Z ro e +a r Je 4> r Q(T)dx (18) 

0 


Y(t) can be approximated using only a subset of the modes for last term in 
eq. (17) and using eqs. (8) and (18), eq. (17) becomes 


( -1 A A-l AT A a -2AT^ 

Y(t) = \K + <D a <D jQ(t) + 1^-K MK + O a <D J Q(t) 

a a -3aT^ aa 

+ \K MK MK + <I> a O jQ(t) + 4>Z(t) (19) 


where the A denotes a reduced set of modes (m < 2n) and Z can be 
calculated by numerically integrating eq. (7) 

This expression can be expanded, assuming a C N differentiable forcing 
function, to 

1 4 


I 


N 

X \(( Y ' 1 -* 1 A A-iAT^iJi-l)] AA 

L^-K mJ K +0 a O J Q(t) J + OZ(t) (20) 

i=l 


Compared to eq. (14), the alternate formulation of the FDM (eq. (20)) does 
not require the solution of a convolution-type integral. In addition, the 
last term of eq. (20) is identical to a mode-displacement solution and as 
such the form of the FDM as given by eq. (20) is more suited for inclusion 
into existing computer codes. 


Second-Order or Natural-Mode Formulation 


Expressions which use the undamped natural modes can be developed in 
an analogous manner and result in an expression similar to eqs. (19) and 
(20). This is accomplished by beginning with the undamped modes of the 
second-order system of equations (eq. (1)). The modes, <j) r , are determined 
by solving the following eigenvalue problem 


K<J> r = <M>r 


( 21 ) 


The modes are normalized as follows 


T 

<j> r M <j) r = 1.0 


T 

<J>r 


K <]) r = (0 


2 

r 


1 5 


so that 



Hence, the displacement response can be represented as 


u(t) = 


X <Mr(t) 


r=l 


( 22 ) 


Using eq. (22), and premultiplying eq. (1) 


by <t> r results in 


2 T 

q + A q+ Q q = [♦] Q(t) 


(23) 


where 

A = [<!> ] T c [<t>] 

and 

tf = [>] T K[<l>] 


where [<{>] is the matrix of undamped eigenmodes, is a diagonal matrix 

2 2 

whose diagonal terms can be represented as Q jj = (0 j and, for proportional 
damping, A is also a diagonal matrix whose diagonal terms can be 
represented as A a = 
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If we assume the forcing function is differentiable, eq. (23) can be 
differentiated twice and back substituted into eq. (23) to produce the 
following expression, which is similar to eq. (16) 


“2 t -2 ”2 r p 

q(t) = Q [(|)] Q(t) - Q AQ [(j)] Q + 


2 -2 -2 -2 - 2 l T " 

fi An AQ - Q Q J(j) Q 


t -2 -2 -2 -2 -2 -2 _2~l (3) 

Q Q A - Q AQ AQ A + Q AQ J q (t) 

r -2 -2 -2 -2 -2] ( 4 ) 

+ [_Q Q - Q AQ AQ J q (t) (24) 


Assuming proportional damping and zero initial conditions, the solution to 
eq. (23) can be written as 


t 

1 f CrC*>r(t-T) T 

q r (t) = Trr J e sin c °dr(t-T)<l)rQ(T) dx (25) 

w drO 


where 

to dr = Vw r 2 - (C r w r ) 2 


If eq. (25) is differentiated four times with respect to t and the expressions 
for q(3) and q(4) are substituted into eq. (24), the entire expression 
reduces to the following 



( A a- 2 aT f *1 '1 A A" 2 * A- 2 aT "\j_ 

= \K - <j) ft <)> jQ(t)-(K CK -<}> ft Aft <|> J ~ 


u(t) = (K - < j > ft <j> 

< -1 -1 -1 -1 -1 
,K MK - K CK CK 


K CK - <j) ft A ft <|> 

). ($ - $ &\s 2 k 6' 2 $ t )|q(,) 


A A 

+ <1> q(t) 


(26) 


Equation (26) agrees with results presented in reference 13 and, as shown 
in reference 13, is also valid for non-proportionally damped structural 
systems. Assuming higher-order piecewise differentiable forcing 
functions, the method above can produce successively higher-order modal 
methods and, as such, is just another formulation of the FDM. The 
expression for an N th -order modal method can be expressed as 





(27) 


where 


B r = 


B I,r‘ 


_ B 2,r. 



-1 -1 
-K CB , r _i - K MB 2 r _! 

B l,r-1 


1 i 

r -n 

>• 

to 

o 

II 

o* 

1 


and 


A r = < i > 


A l.r 
- A 2,r. 


aT a 

< i > = 4 > 


-ft Aa i - ft 2 A 2>r .i 

A 1 ,r- 1 


<■2 


aT a 

<j> , A 0 = 4» 


ft 

L 0 J 


aT 

I* 
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The Dynamic-Correction Mmethod (DCM) 


The dynamic-correction method (DCM), reference 13, assumes a solution to 
eq. (2) in the form 


Y(t) = Yp(t) + Y c (t) 


(28) 


where Yp(t) is a particular solution of eq. (2) and Y^(t) is the 
complimentary solution which represents the effects of initial conditions. 
In modal form, Yp(t) and Y^(t) can be represented as 

Y P (t) = [<t>] Zp(t) 

and (29) 


YcW = [«>] ZcW 

where Zp(t) and Z^(t) are the vectors of particular and complimentary 
solutions to the modal coordinate equations (eq. (7)). 

Using eqs. (6) and (28) 


Y(t) = [<t>] Z(t) = [<K][Z(t) - Zp(t)] + Yp(t) (30) 


The fundamental principle of the DCM is that if you have an exact 
particular solution to eqs. (2) and (7), you can approximate the response 
(eq. (30)) using a reduced set of modes as shown below 


aa r a a "1 
Y(t) = <DZ(t)+[Yp(t)-<I>Zp(t)J 


( 31 ) 
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It can also be shown that in the limit as N goes to infinity, two terms in 
equation (20) can be written as 


lim 

N-*oo 


A 

-<x> 


' N 

X 

Lr=l 


A*r 

ot O Q(t) 


A A 

= OZp(t) 


and 


(32) 



r 

‘ N 

> 

< 



> 

lim = 


_£,U-K M K * Q(t) 


N-»oo 

< 

L r=l J 



= Y P (t) 


Hence, if an infinite number of integrations-by-parts are assumed in the 
FDM or if the convolution integral vanishes (e.g., for a polynomial forcing 
function of a lower order than the order of the FDM) the FDM would be 
equivalent to the DCM of reference 13. Also, if an exact solution to the 
convolution integrals of eqs. (8), (12), or (14) exists and is used; the 
response can be calculated without errors caused by the approximation of 
the forcing function. For the case of a sinusoidal forcing function, 
equations (32) will converge as shown above provided the frequency of 
the highest mode used in the approximation is larger than the forcing 
function frequency. 
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RESULTS AND DISCUSSION 


Structural Example: Two-Degree-Of-Freedom Problem 

A simple, two-degree-of-freedom spring-mass problem (fig. 1) with a 
sinusoidal forcing function was analyzed to compare the accuracy of the 
MDM, MAM, FDM, and DCM. This problem was also investigated in 
reference 13 and included in that reference are the particular solutions for 
polynomial as well as sinusoidal forcing functions. As shown in figure 1, 
the sinusoidal forcing function, sin((Oft), is applied to the second mass. The 
natural frequencies are (0j= 19.54 rad/s and © 2 = 51.17 rad/s. The system 
is proportionally damped (A (eq. (23)) is diagonal) if a = K 1 /K 2 (ref. 15). For 
the stiffnesses chosen, this corresponds to a value of a = 1.0. The accuracy 
of each method is assessed by a time-integrated error norm, which is 
defined as 


J I u i(t) - u*j(t)| d t 

£j (%) = X 100 (33) 

J I Ui(t) I dt 

0 


where uj(t) is the calculated response using all the modes and uj(t) is the 
approximate response using a subset of the modes. Results were 
calculated using both the real and damped modes (eqs. (20) and (27)) 
respectively. For this problem, the FDM used was of order four (N=4 in eqs. 
(20) and (27)). The time, X, selected for integrating the error was chosen to 
be X = 16Jt/(Of. 

Results of the error as a function of the forcing frequency for the 
undamped modal solution using one real mode for the proportionally- 
damped case (a =1.0) is shown in figure 2. In general, the forcing function 
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frequency must be lower than the highest natural frequency used in the 
approximate modal response for accurate results. As shown in figure 2, 
the accuracy increases as the order of the modal method increases. The 
results for the lower range of frequencies are shown more clearly in 
figure 3, which is an expanded error scale of figure 2, that the FDM (N=4) 
and DCM are similar and more accurate than the lower order methods such 
as the MDM (N=0) and the MAM (N=l) for (Of < 20 rad/s. For cof > 20 rad/s, 
the DCM remains slightly more accurate than the FDM, however, as the 
forcing frequency approaches the second natural frequency ((0f= 51.17 
rad/s) all methods produce inaccurate results. It should be noted that 
results using one real mode or two damped modes are identical for the 
proportionally-damped case. 

Results for the non-proportionally-damped case (a = 20) using two 
damped modes, are shown in figure 4. Results are similar to the 
proportionally-damped case with the exception that the DCM exhibits 
surprisingly good results at forcing function frequencies close to the second 
natural frequency. This result is unexplained at present and is believed to 
be fortuitous and, hence, it is recommended that all modal methods should 
include modes whose frequencies bound the frequency of the forcing 
function. A comparison of the damped-mode solution ( using two damped 
modes (eq. (20)) and the undamped solution (using one real mode (eq. 

(27)) is shown in figure 5 for the non-proportionally-damped case (a = 

20). As shown in figure 5, the damped-mode solution using two damped 
modes (dashed lines) produces more accurate results than those using only 
one real mode (solid lines). The damped-mode solution for the FDM and 
DCM are nearly equivalent and result in the smallest error for frequencies 
as large as 30 rad/s. Hence, it may be beneficial, in some cases, to use the 
damped modes to obtain a more accurate solution. 

The FDM produced similar results to the DCM for forcing frequencies below 
the first natural frequency. A comparison of the modal methods for a 
forcing frequency (Of = 10 rad/s is shown in figure 6. Once again, the 
higher-order modal methods result in more accurate solutions. The large 
relative errors near X = 0 are due to the zero initial conditions which cause 
the denominator of eq. (33) to approach zero at X = 0. As explained in 
reference 12, the increase in accuracy with the order of the modal method 
is due to the addition of terms which are functions of the generalized 
stiffness and mass matrices and the force vector and its time derivatives. 
These additional terms approximate the effect of the higher modes which 
were neglected in the modal summation. 
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Thermal Example: Rod Heated At One End 


Results of references 10 and 11 indicate that Lanczos vectors can be 
effective reduced basis vectors for solving linear and non-linear thermal 
problems. Since the accuracy of the Lanczos vectors is comparable to that 
of the MAM for structural dynamic problems, it was expected that higher- 
order methods, such as the MAM, FDM and DCM, would be effective in 
solving complex thermal problems. The problem selected to study is 
similar to that presented in reference 10 with one exception; the present 
problem assumes the temperature at the right end of the rod is 
constrained to zero. The forcing function is a ramp up from zero to a peak 
value at time t = 10 sec and a ramp down to zero at time t = 20 sec as 
shown in figure 7. The error function used to evaluate convergence is a 
spatial error norm similar to that used in reference 12, namely 


V ( T - T a 7 ( T - T a ) 

^ 1 ( 34 ) 

A total of twenty equally-spaced finite elements were used to model the 
problem. The temperature distributions in the rod as a function of the 
number of modes for the MDM, MAM, and FDM are shown in figures 8a to 
8c for time t = 10 sec. Because the forcing function is linear, the FDM 
(having an order of 2) is exactly equivalent to the DCM. The FDM 
converges with 5 modes to the exact solution as compared to 8 modes for 
the MAM and 18 modes for the MDM. A spatial error norm similar to that 
shown in eq. (34) is used to compare each method for time, t = 10 sec (fig. 
9). The effectiveness of using higher-order modal methods for reducing 
the size and computational effort of thermal problems is illustrated in 
figure 9. The FDM or DCM require about 28-percent of the number of 
modes as compared to the MDM and about 63-percent of the number of 
modes as the MAM for an accurate thermal response. 



CONCLUDING REMARKS 


The present study extends the development of higher-order modal 
methods to include non-proportionally damped systems using both first- 
order (damped modes) and second-order (natural modes) modal- 
superposition methods. The study compares the accuracy and convergence 
of the mode-displacement, mode-acceleration, force-derivative, and the 
dynamic-correction methods (MDM, MAM, FDM, and DCM respectively) in 
solving proportionally and non-proportionally damped structural dynamic 
problems as well as transient thermal problems. The higher-order modal 
methods, such as the FDM or DCM, are very effective in solving structural 
dynamic problems. The damped mode solutions are found to be effective 
in solving a non-proportionally damped two-degree-of-freedom problem. 

Results of a two-degree-of-freedom spring-mass-damper system indicate 
that, for the proportionally damped problem, a solution using two damped 
modes produces identical results as one using only one natural mode. 
Hence, there is no advantage in using a damped modal solution to solve a 
proportionally damped problem. However, for the non-proportionally 
damped problem, the use of two damped modes produces more accurate 
results than the single natural mode case. The DCM has the lowest 
percentage error of all the mode-superposition methods over the 
frequency range of 2 to 50 rad/s. The FDM, having an order of four (four 
integrations-by-parts (N=4)), produced similar results to the DCM up to a 
forcing frequency of about 35-40 rad/s. For the proportionally damped 
problem, all the methods were inaccurate near a frequency of 50 rad/s 
(close to the second natural frequency of the system). 

A one-dimensional heat conduction problem (rod heated at one end) was 
also investigated to evaluate the usefulness of using the FDM to solve 
thermal problems. The higher-order modal methods such as the FDM were 
very effective in solving the thermal problem. Until now, modal methods 
were inefficient in solving thermal problems because the nature of the 
problem required the inclusion of the higher modes for an accurate 
solution. The ability of the FDM to approximate the effects of the higher, 
but neglected, modes resulted in an accurate solution using only five 
modes out of a total of twenty as compared to the MDM which required 18 
modes for an accurate solution. 
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Two-degree-of-freedom spring-mass problem. 



Modal Methods (1 real mode 
or 2 damped modes) 

Mode displacement (MDM) 
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Figure 2.- Comparison of modal methods as a function of forcing frequency for 
a proportionally-damped two-degree-of-freedom problem (a = 1). 
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Figure 3.- Comparison of modal methods as a function of forcing frequency for 
a proportionally-damped two-degree-of-freedom problem (a = 1). 



Modal Methods (2 damped modes) 



Ui 


30 


Figure 4.- Comparison of. modal methods as a function of forcing frequency ^ 
non-proportionally-damped two-degree-of-freedom problem (a - zu) 
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Figure 5.- Comparison of damped and natural mode solutions for a non- 

proportionally-damped two-degree-of-freedom problem (a = 20). 




x (sec) 

Figure 6.- Comparison of modal methods as a function of time for 
proportionally-damped problem (a = 20) with (Of -10. 
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One-dimensional haet conduction problem: rod heated at one end. 
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Figure 8.- Temperature distribution along a rod heated at one end at t 10 sec 



Time = 10 sec 

Mode-acceleration method (MAM) 



b) Mode-acceleration method (MAM) 





TEMPERATURE DISTRIBUTION ALONG ROD 

HEATED AT ONE END 
Time = 10 sec 

Force-derviative method (FDM) 
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c) Force-derivative method (FDM) 



= (T-T a ) T (T-T a ) 
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Figure 9.- Convergence of modal methods for a one-dimensional, transient 
thermal problem at t = 10 sec. 
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